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Abstract 

In this paper, we constrain CP violation in the Higgs sector using the measured 
signal strengths in the various Higgs search channels. To this end, we introduce a 
general parameterization for a resonance which is an admixture of a CP-even Higgs- 
like state and a CP-odd scalar. By performing a fit to the available data from the 
Tevatron and LHC experiments one obtains constraints on the mixing angle and the 
couplings of the resonance to Standard Model fields. Depending on the couplings, 
sizable mixing angles are still compatible with the data, but small mixing is in general 
preferred by the fit. In particular we find that a pure CP-odd state is disfavored by the 
current data at the 3<r level. Additionally we consider a mixed fermiophobic resonance 
and a model with two degenerate mixed resonances and find that both scenarios can 
successfully fit the data within current errors. Finally, we estimate that the mixing 
angle can be constrained to a < 1.1 (0.7) in the full 8 TeV (14 TeV) run of the LHC. 



1 Introduction 



Recently, the ATLAS and CMS collaborations have reported on the discovery of a new 
bosonic resonance with mass in the range 125-126 GeV at the Large Hadron Collider (LHC) 
[TJ[2], which has been corroborated by an excess observed by the CDF and D0 experiments 
at the Tevatron pUE]. While the current data is in agreement with expectations for the 
Standard Model (SM) Higgs boson, the experimental uncertainties are still large, and thus 
other possibilities still need to be considered. 

Of particular interest are the spin and CP quantum numbers of the new particle. Since 
it is known to decay into photon pairs, it cannot be a spin-1 particle. A spin-2 resonance 
may be distinguished from a spin-0 resonance by analyzing angular distributions in the 
77 0E], ZZ* U [M], WW* Ivtv 0, and Z 7 U-f PHU] decay channels, or angular 
correlations in associated production with jets Furthermore, a spin-0 particle <fi may 
be CP-even, CP-odd, or a general mixed CP state. The CP properties can be determined 
from angular distributions in ZZ* — > 4£ [7JE], angular distributions of the jets in + 2 jets 
production [12], or spin correlations in (f> — > t + t~ decays [13] . 

However, the analysis of distributions becomes viable only if a sufficient number of events 
has been accumulated in a given channel. At this early stage, however, one can already 
constrain the CP properties from the observed production rates and decay branching frac- 
tions p3UT5] . This mainly follows from the fact that a CP-odd even-spin particle cannot have 
renormalizable tree-level couplings to two gauge bosons. Based on this approach, Ref. [15] 
finds that a CP-odd pseudoscalar is disfavored compared to a CP-even scalar, although their 
conclusion is not based on a global fit to the known data and thus difficult to interpret 
statistically. The goal of this paper is to carry out such a fit, in a general setup where the 
125-GeV resonance can be an arbitrary mixture of CP-even and CP-odd components and 
can have modified couplings to SM fermions, as well as new couplings to SM gauge boson 
mediated through higher-dimensional operators. 

The model setup is explained in more detail in section [2j The possibility of general CP 
mixing leads to modified decay branching fractions and production rates, which is discussed 
in section |3j In section HI these observables are then confronted with the available experi- 
mental data from July 2012 to put constraints on the amount of CP mixing and coupling 
parameters. Finally projections for how these bounds may improve with additional data 
from the LHC are presented in section El before concluding in section EJ 

2 Setup 

Throughout this paper, it will be assumed that the 125-GeV resonance observed by ATLAS 
and CMS, denoted 0, is a scalar, but its CP properties are left unconstrained. In general, it 
can be a mixture of a CP-even Higgs-like scalar H and a CP-odd scalar A: 

<fi = cos a H + sin a A. (1) 

CP mixing in the Higgs sector can appear in many extensions of the SM. Two of the simplest 
possibilities are a complex singlet extension of the SM [16] or the Two-Higgs-Doublet Model 
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(THDM) pig. As a result of eq. (GQ), the tree- level couplings of to W and Z bosons 
are reduced by a factor cos a compared to the SM, since AW + W~ and AZZ couplings can 
be realized only through operators of dimension five or higher and thus are expected to be 
generated through loops of heavy particles. 

Generically, we assume that the orthogonal state 

(f)' = — sin a H + cos a A (2) 

is much heavier than cf) and evades current search limits through its modified couplings 
compared to the SM Higgs. 

In a general extension of the SM, the Yukawa couplings of the CP-even and CP-odd 
components of (f) are free parameters. However, existing data on the fermion masses and 
mixings essentially demands that the up-type and down-type Yukawa matrices can be written 
as the SM Yukawa matrices Y u,d ' e times some overall constant for each matrix. This is 
described by the Lagrangian 

£vuk = - VuY^UiUjH - ydY^didjH - y&f-fiitjH 

— ix u Y^UiUjA — ix^Y^didjA — ix^Y^IiijA + h.c. , 

where y U)( j and x Uj( j parametrize the strength of the CP-even and CP-odd Yukawa couplings, 
respectively, relative to the SM coupling strength. In particular, the THDM types I and 
II fit in this pattern. Note that the framework in eqs. ([1]) and (J3]) is general enough to 
accommodate the possibility that H itself is a mixture of several CP-even states — in this 
case a, y u and would be functions of the 3x3 Higgs mixing matrix. The SM corresponds 
to the choices a = 0, y Ui d = 1, and a; U) d = 0. 

While there are no renormalizable couplings of the CP-odd component A to the SM gauge 
bosons, higher-dimensional interaction operators may be induced through loop corrections 
of heavy new fields. In an effective field theory formulation these interactions are given by 

£ dim5 = \ 1 ^ r AG IBf Cr + \-^—AB, v B^ + \^—AW v Jw^ , (4) 

with G^ u = e^ ua ^G a j3, etc. and v = 174 GeV is the electroweak vacuum expectation value. 
The normalization is chosen such that the Feynman rules have a prefactor Ci/(16ir 2 v), see 
appendix. We assume that the coefficients originate from perturbative new physics, so that 
q < 4tt. We do not consider dimension-five operators for the coupling of the CP-even 
component H to gauge bosons, since the effects are typically small compared to the tree- 
level HWW and HZZ couplings, while the loop-induced if 77 and H gg interactions can be 
sufficiently generally described by the modified Yukawa couplings in eq. (j3J). 

*The LHC Higgs data has been analyzed in the context of specific realizations of the THDM in several 
recent papers [TS] . 
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3 Decay Widths and Production Rates 

Let us begin by disregarding the dimension-5 operators in eq. (pE}, to illustrate the effect of 
the CP mixing and modified Yukawa couplings. Compared to the SM, the partial widths for 
the tree-level decays are given by 

T[<f>->WW*] Tty^ZZ*] 2 



T SM [H^WW*} T SM [H^ZZ*] 



cos a , (5) 



r [(j) T+T 2 2 

- — — — — T = {y d cos a) + {x d sm a) , (6) 

T SM [H T+T \ 

r[(j) -> cc] = cog 2 + RCB gin 2 ^ ^ 
r SM [-H -> ccj 
ru -a mi _ . E 

(ydcosa) 2 + bb (a; d sina)% (8) 



r SM [# 66] 

while the loop-induced decays widths read 

99] 2 \yuH 1/2 (T t ) + y d H 1/2 {T h )\ 2 . 2 \x n A l/2 (n) + x d A l/2 {T h )\ : 
— cos a — — - — -. — -r^r h sm a K yy - 
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(9) 



T SM [H^gg] \H 1/2 (r t )\ 2 \H 1/2 (r t ) 

r[0 -> 77] 2 l|yu-gi/2(n) + ggagi/gfa) + VdH l/2 {T T ) - #i(t w )| 2 , . 

r SM [^^77] C ° S " |f^ 1/2 (r t )-i/i(rw)| 2 1 j 



+ sm a 



3^1/2 

, ||x u Ai /2 (r t ) + |x d A 1/2 (r b ) + x d A 1/2 (T T )[ 



where Tf = m^/(4m^), and 



M M ( r - 1 )^ r ) + r MM 3(2r-l)/(r) + 3r + 2r 2 /(r) 
H 1/2 (t) = — , H 1 [t) = — 2 , A 1/2 {t) = — — , (11) 

{arcsin 2 (- v /r) (r < 1) , 

x a i +v ^i77 . y (12) 

For the SM decay rates into gg and 77 in (jUJ) and f lTUj) . the contributions from tau leptons 
and bottom quarks may be safely neglected, but for decays they can be enhanced by large 
Yukawa factors y d and x d and thus need to be included. The factors R x incorporate the dif- 
ference between the QCD corrections for scalar and pseudoscalar decays, R x = 1^™^^] ■ 
They deviate from unity by only small amounts of less than 1%. 

The production rates at the Tevatron and LHC for final states X = WW* , Z Z* , 77, rr 
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can then be written as 



<?sm\pV ->■ H X] 



SM ,, rl ( ' ;,) ) 



Jgs^ T77 — — ~T + JVBF^ TFT — . 777777Ii ) x Ti x 



T SM [F^^] JVDP T SM [H^WW*]J T^ tot r SM {H^X] 

Here / 99 and /vbf denote the fractions with which the two dominant production modes, gluon 
fusion and vector-boson fusion, contribute to the total production cross section, respectively 
(fgg + /vbf = 1)- T^tot is the total decay width of 0, which is given by 

r ,tot ~ r[0 -> ww*\ + r[0 -> zz*] + r[0 -> 66] + r[</> -> cc] + r[0 ->• r+r-] 

to very good approximation^ (similar for rf^[ ot ). For the SM partial widths we take the 
values from Ref. [3]. 

For the bb final state, the two leading production modes are not experimentally viable 
due to large backgrounds. Instead, the experimental collaborations focus on associated 
production with a W or Z gauge boson, which scales with cos 2 a according to eq. (jSJ). Thus 

2 rf, M tot r[0 -» 66] 

r 6 g = cos a x — ! — x . (15) 

r^tot r SM [-ff-^66] 

If the contributions from new-physics induced higher-dimensional operators in eq. (jl]) are 
sizable they lead to additional contributions to the partial widths into gauge-boson pairs. 
The relevant Feynman rules are listed in the appendix. Including these terms, one obtains 
(with the approximation R" m 1) 

r[0->##] 1 , , | , , 

cos a y u -Hi/ 2 (T t ) + yd-ni/2(Tb) 



T SM [H^gg) \H 1/2 (r t )\ 2 \ |OT w " n i/JV D/| (16) 

+ sin 2 a |a; u A 1/ 2(rt) + x d A 1/2 (r b ) + \Z2c G /# 2 | 2 J , 



r[0^ 77 ] 1 



r SM [# -> 77] ||%(T-t) - #i(r w )| 2 

x ^cos 2 a | |y UJ ffi/ 2 (rt) + §2/d#i/2(Tb) + y&Hx/ifa) - #i(t w ) | 2 (17) 

^u^l/2(n) + |^d^l/2 (^b) + ^dA/2(T r ) + ' 



+ sin 2 a 



V2e 2 



where cq == cos^w and sg == sin# w are the cosine and sine of the Weinberg angle, 
respectively. For the jZ decay one arrives at a similar expression, which we do not write 
down here since it is rather lengthy [19] - In fact, the SM contribution to this decay channel 



^We do not consider possible non-standard decay channels of <fi in this paper. 
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is rather small and thus irrelevant for the current early stage of Higgs searches. However, the 
dimension-5 operators in eq. fll]) could potentially lead to a much larger result that would 
dominate over the SM contribution, in which case one can write 



2 {„___ „_\2 f m 2 _ m |^3 



r[0 ->■ 7Z] as sin 2 a 



?(47r) 5 f 2 m; 



For the four-body decay modes mediated by WW and ZZ pairs one finds, using CalcHEP 
for the numerical phase-space integration, 

T\d> -> WW*1 9 9 cL 

cos 2 a + sin 2 a— ^- x 0.155, (19) 



t sm [h -> ww*] (4tt; 

2 2 (s^Cs + dcH/) 2 , . 

- cos 2 a + sin 2 a 9 * x 0.074 . (20) 



r SM [if->ZZ*] (4vr) 

Finally it is important to note that there are no interference effects between the CP-even 
and CP-odd contributions in the inclusive rates, in contrast to specifically CP-sensitive 
observables, such as certain angular distributions [7]l8tlT2|fT3]. 



4 Numerical Analysis of Summer 2012 Data 

In this section, the formalism of the previous two sections is applied to the experimental Higgs 
search results released in July 2012 by the Tevatron and LHC collaborations p~M3 j. l21| . l22j . It is 
shown that the CP properties of the new resonance can already be constrained substantially 
from our current knowledge of its production rates and branching fractions, even though the 
experimental uncertainties are still large. The relevant input data, as read off from the plots 
in Refs. pH3l l2T|[22] . is summarized in Tab. [U The following values are taken for the relative 
production rates: 

• Tevatron inclusive: f gg « 0.78, /vbf ~ 0.22 [3]; 

• LHC inclusive: f gg w 0.9, / V bf ~ 0.1 [21]; 

• LHC 77 VBF enhanced: f gg w 0.25, / V bf « 0.75 [21]; 

• LHC tt analysis: f gg ~ 0.5, /vbf ~ 0.5; 

estimated from the observation that inclusive and VBF-enhanced measurements of this 
channel contribute with approximately similar significance [23] . 

Here "VBF enhanced" refers to Higgs searches with a set of cuts that enhance the relative 
contribution of the vector-boson fusion production mode, characterized by two energetic jets 
with a large rapidity gap. Since both these searches and the inclusive LHC measurements 
receive contributions from gluon fusion and VBF, there is some degree of correlation between 
them, which is taken into account with a covariance matrix in the x 2 fit- 
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Experiment 


X 


rx/r SM 


Ref. 




WW* 


1.24 ±0.45 






zz* 


1.39 ±0.60 






bb 




m 


ATLAS 


TT 


0.45^ 






77 


1.79 ±0.50 






77 (VBF enh. 7 TeV) 


4.19 ±2.10 


EI] 




77 (VBF enh. 8 TeV) 


1.24 ± 1.57 






WW* 


RQ+ - 46 






ZZ* 


n 79+0.48 
u - ' z -0.35 




CMS 


bb 


0.48+°; 8 72 


m 


TT 


0.08i°;?i 






77 


1.56 ±0.47 






77 (VBF enh.) 


2.30 ± 1.26 


[22] 




WW* 


32 +113 

U - OZ -0.32 




CDF/D0 


bb 


1 O7+0-74 
±>y ' -0.68 






77 


o £9+2.96 
°- uz -2.54 





Table 1: Experimental results for Higgs production rates in different final-state channels 
from Tevatron and LHC used in this analysis. Separately shown are the values for 77 final 
states with cuts to enhance the VBF production mode from ATLAS and CMS. 

Fitting this data to the SM predictions we find 

X 2 SM = 13.3 (16d.o.f.). (21) 

For a x 2 distribution with 16 degrees of freedom, the 68% (95%) confidence limit (C.L.) 
corresponds to \ 2 = 17.0 (25.0). Thus one can see that the overall agreement of the data 
with the SM prediction is very good. For the two parameter plots that will be shown later 
the 68%, 95% and 99.7% contours correspond to x 2 = 15.9, 23.7, and 32.9. 

4.1 Single resonance 

We first consider the scenario described by equations ([T]) and OH])- Here both H and A have 
renormalizable couplings to the SM fermions, while the couplings of A to SM gauge bosons 
are induced at the one loop level through SM fermion loops. We assume that there are no 
other new physics states that generate couplings of A to SM gauge bosons, and therefore 
cg — c b = cw — in 01]). The fermion couplings of H and A are allowed to deviate from 
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Figure 1: Quality of the fit to experimental 
data in the x u — a plane, for x^ = and 
Vu — Ud — 1- The orange (grey) shaded areas 
agree with the data at the la (dark), 2a 
(medium) and 3a (light) level. The blue star 
shows the best fit point, while the black dot 
corresponds to the SM. The solid (dashed) 
lines are contours of constant r 77 {r Z z)- 



their SM values through the parameters y Ut d and x U) d respectively. As mentioned in section [2, 
it is assumed that = 125 GeV <C m^. 

Let us first consider the case where the CP-even state has SM-like Yukawa couplings, 
Vu — Vd — 1) while the CP-odd couplings x Ut d are unknown. In the limit of zero CP-odd 
couplings, all channels are uniformly suppressed by cos 2 a. For this particular point, we find 

a < 0.76 (95% C.L.) , (22) 

while the best fit is for a = 0. When allowing the CP-odd couplings to float freely, the overall 
rate suppression from the mixing can now be offset with an increase in the production rate 
when x u > 0. In fact, large values of x u are favored in the fit, so to ensure perturbativity 
of the top Yukawa coupling, we impose an upper limit x u < 3. The effects of x& are more 
subtle. It can increase the total width and thus suppress all but the rr and bb channels, so 
large values of Xd do not produce a good fit. Overall, we find that a marginally better y 2 
than for the SM is obtained for nonzero but small mixing a = 0.07, maximal x u = 3 and 
vanishing x&. 

The large value of x n together with a small mixing leads to slightly enhanced signal rates 
at the 10% level across all channels, which is slightly favored by the current data. The overall 
quality of the fit in the x u -a plane for Xd = is shown in Fig. [T]on the left. Mixing angles 
of up to a = 1.3 are compatible with the data at the 95% C.L. level. 

Close the a = %/2 the field <j) becomes mostly CP-odd, and the signal rates r Z z and 
r ww become strongly suppressed. For smaller mixing, both r zz and be enhanced 

or reduced relative to the SM, however an enhancement of the di-photon rate by more than 
50% is only possible outside of the la region. 

Let us now consider the case where the CP-even Yukawa couplings can vary with respect 
to the SM. Due to the additional free parameters, the predicted rates in the different channels 
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are less strongly correlated with each other. For nonzero mixing, there is some redundancy 
in the couplings x u , y u and x^, yd that lead to almost flat \ 2 distributions in some directions. 
Numerically we find the minimum of the \ 2 distribution at 

a = 0.16, a; u = 3.0, x d = y n = , y d = 0.74. (23) 

With x 2 = 7.6 this parameter point lies significantly below the SM fit, however at the 
expense of having five additional model parameters. Both the CP-mixing scenario and the 
SM have x 2 /(d.o.f.) < 1, i.e. the current data does not conclusively favor one model over 
the other. Nevertheless we will indicate the most preferred regions with A% 2 < 1 from the 
minimum in the plots. Should future more precise data have similar central values, those 
regions would be strongly favored. 

The good quality of the fit can be understood as follows. First, the gluon fusion produc- 
tion is reduced to roughly 65% compared to the SM, which is mostly due to the vanishing 
CP-even Yukawa coupling, y u = 0, while the small mixing angle suppresses the coupling of 
the CP-odd component to gluons. The mixing also slightly suppresses the VBF channel to 
about 85%. The di-photon decay width is naturally enhanced since the destructive inter- 
ference between the top and the W^-boson loop goes away. Additional enhancement comes 
from the large decay width of the CP-odd component into photons, and from a reduction 
of the dominant bb width due to yd < 1. The change in the WW and ZZ channels is more 
balanced, since the increased branching ratio to those final states is compensated by the 
overall reduced production cross section. Finally r TT is suppressed by about 50% due to the 
smallness of yd- 

The dependence of the fit on the different model parameters is illustrated in Fig. [2j The 
various plots show a scan over two parameters while the remaining free parameters are set to 
the best-fit values in f l2"3"j) . We can easily see that the best-fit regions typically have r 77 « 1.5 
while the ZZ and WW rates are kept closer to one. 

The correlation between x u and the mixing angle is very strong. When x n is reduced a 
larger mixing is required to fit the data, since otherwise the total production cross section 
becomes too small. Similar to Fig. [1] the ZZ and WW channels are strongly suppressed for 
a > 1.0, such that this region never leads to a satisfactory fit. 

The interplay between y d and a in the top right plot of Fig. [2] is again more subtle. 
Since x n is fixed and y u = here, increasing yd leads to a suppression of the gauge boson 
channels as the total width goes up. To some extent this can be compensated with an 
increase in a which increases the total production cross section. Eventually this leads to a 
strong enhancement of r TT such that regions above yd ~ 2.5 are excluded here. 

The last two plots illustrate the redundancy in the couplings x u -y u and £d~2/d- In the 
bottom left plot of Fig. [2j one can see that within the la contour one can trade x u for y u , 
with the ratio of the two couplings roughly given by the mixing angle. The preference for 
smaller y n comes mostly from r 77 , since y n < 1 reduces the destructive interference and thus 
increases the decay rate of the CP-even component into photon pairs. Thus the redundancy 
between x n and y u can eventually be broken with more precise data on r 77 . The quarter- 
circle that is described by the fit contours in the Xd — yd plane can easily be understood 
from ([8]): it corresponds to contours of constant bb width. These couplings are only weakly 



8 



3.0F 



2.5 



2.0 



1.5 



1.0 



0.5 



0.0 




0.0 



0.5 1.0 

a 



1.5 




1.5 



yy 




2.0- 








1 

1 




1 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 



2 4 6 8 10 

Xd 



Figure 2: Quality of the fit for the single resonance model, in the a-x n (top left), a^tjd (top- 
right), x u -y u (bottom left) and £d^/d (bottom right) plane. The remaining free parameters 
are set to their best fit values, see text for details. Colors and contours as in Fig. d In 
addition the blue (very dark) shaded region indicates Ax 2 < 1 relative to the best fit, and 
the dotted lines in the top-right plot are contours of constant r TT . 



9 



constrained from the ^77 and 4>gg couplings due to the smallness of the bottom Yukawa 
coupling, so that this redundancy is difficult to resolve in general. 

4.2 Two near-degenerate resonances 

As pointed out above, large CP mixing can be realized in models with several scalar multi- 
plets. However, in the context of concrete models there are tight constraints on CP mixing 
in the Higgs sector from electric dipole moments [24J. It was shown in Ref. [21] that these 
bounds are substantially relaxed for a near-degenerate Higgs spectrum, i. e. if the two or- 
thogonal states and 0' have almost equal masses m^. 

If |m^— 7710/ 1 < 1 GeV, both states would contribute to the resonance observed by ATLAS 
and CMS. This scenario is explored in more detail in this subsection]*] The relevant branching 
fractions and production rates for 0' can be derived from the formulae in section [3]by making 
the appropriate replacements of the mixing angles. The observed rates are then given by 

_ a y-^ X ] + *\p$ -> ^ -> x) 

TX a m \p¥ ->H -> X] ' 1 } 

The signal rates for 0' are obtained from the rates by a shift a — > n/2 — a. It therefore 
follows that the combined rates rx are symmetric under this transformation. For definiteness, 
we choose a < tt/A when searching for the best-fit point. 

Letting all model parameters float freely, we find the minimum of the \ 2 distribution at 

a = 0.38, x d = y u = 0, x u = 0.57 , y d = 0.75 . (25) 

The quality of the fit is marginally better than in the single resonance model. In both 
scenarios, the di-photon channel is enhanced for y u = and x u > 0, due to absence of 
destructive interference as explained in section 14.11 resulting in r 77 1.6 for the best- 
fit point. In addition to the enhancement of r 77 as before, now there is also a relative 
enhancement of the VBF channel compared to the gluon fusion production, because of the 
contributions of both and 0' exchange, such that we obtain r^y BF > 2. Furthermore, 
there is a slight suppression of rzz and rww from the mixing, and a stronger suppression 
of r TT due to y d < 1. Altogether this leads to a very good fit to the data. Comparing with 
the single resonance model, the preferred value for x u is now much smaller and the mixing 
angle is increased, while x d and y Ujd are roughly the same. 

The parameter dependence of the fit in the two-resonance model is illustrated in Fig. |3J 
Compared to the model with a single mixed resonance, there are some marked differences. 
First we note from the left plot that the degeneracy between a and x u is now broken, and 
x u < 1.2 at 95% C.L. Interestingly, this constraint is mainly due to a too strong enhancement 
of r 77 for larger values of x u . On the other hand, a is essentially unconstrained now, except 
that very small mixing angles are disfavored. The latter follows from y u = 0, which leads to 
a strong suppression of rzz and rww when the mixing becomes too small. 

*A similar study, but for Higgs mixing with a CP-even resonance, can be found in Ref. [25) . 
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Figure 3: Quality of the fit for the double-resonance model, in the a-x n (left) and x u -ya 
(right) plane. The remaining free parameters are set to their best fit values, see text for 
details. Colors and contours as in Fig. [I] and 

The t/d dependence shows an inverted behavior. Values above 1.6 increase the total width 
and thus suppress both rzz and r 77 , while very small values lead to a too strong enhancement 
of both channels. Together this suggests some correlation in the x u -yd plane, as can be seen 
in the right plot of Fig. [3J The region that leads to good agreement with the data is relatively 
constrained. The modifications of r 77 and rzz cancel along certain diagonal directions, but 
they start deviating from each other for larger values of the couplings. In addition r TT grows 
with x u and, to a lesser extent, with yd, such that very large values of both couplings are 
also disfavored by this observable. 

The correlations in the Xd-yd and x u -y u planes are very similar to the case of the single 
resonance model, so we do not show them separately. Overall we find that the double- 
resonance model imposes stronger constraints on the variations of the Yukawa couplings, 
while the mixing angle is less constrained in this scenario. In addition the possibility to 
enhance r 77 without modifying the Yukawa couplings of the CP-even component sets this 
model apart from the single-resonance case. 

4.3 Effective theory including new dimension-five operators 

Now we turn to the discussion of the effects of higher-dimensional operators which couple 
the CP-odd component A to gauge bosons, see eq. (BJ, and which are induced by loops 
of heavy new particles. The contribution of these operators, with coefficients cg, cw, and 
Cb, comes in addition to contributions from top and bottom quark loops, which can induce 
sizable couplings between A and 77 and gg pairs through the couplings x Uj d. [The effect of 
top and bottom loops is generally negligible for ZZ and WW pairs, which are dominated 
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Figure 4: Quality of the fit in the effective operator model, in the c^-cb (left) and c-w-oi 
(right) plane. The remaining free parameters are set to their best fit values, see text for 
details. As in the previous figures, the shaded areas agree with the data at the level of 3a 
(light), 2a (medium), la (dark), and Ay 2 < 1 (blue/very dark). The solid lines are contours 
of constant r 77 , while the dashed lines denote constant r^ 7 (left) and rzz (right). 

by the tree-level HZZ and HWW interactions.] 

The expressions for the partial decay widths including the contributions from both the 
modified top/bottom loops and cg,w,b are shown in eqs. (FIdT) - (|211 . To avoid redundancy 
with the results of the previous sections, we set in the following x u = Xd = and y u = Ud = 1- 
In the limit a = n/2 we recover the model [15], a pure CP-odd singlet scalar that only couples 
to the SM through dimension five operators. 

Letting the mixing angle and the coefficients of the dimension five operators float, we 
find the minimum of the v 2 distribution at 

a = 0.28, c G = 0, cb = 0.32, c w = ll.b. (26) 

With x 2 — 8.86 this provides a better fit than the SM but slightly worse than the scenarios 
considered in Sec. 14. H and 14.21 The dominant effect is an enhancement of the di-photon rate 
to r 77 ~ 1.7 from the dimension five operators. The small mixing angle together with cq = 
suppresses all other rates by roughly 10% compared to the SM. 

Since only r 77 is notably modified, the couplings cb and cw are not probed separately but 
through the combination cb + cot 8 2 cw- This is immediately obvious from Fig. H] where the 
iso-contours of r 77 in the cb~cw plane are straight lines. To probe the couplings individually 
one would have to measure the ratio r Zl , which is proportional to cb~cw A combined 
measurement of both rates would then single out a circular region in the cg-c^y plane, 
however the required precision in r^ 7 can only be achieved with higher luminosity at the 
14 TeV LHC [5B]. 
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In the absence of a result for r Zl we can set cb = without loss of generality, and 
analyze the constraints on the mixing angle in this model. In the right plot of Fig. H] it can 
be seen that very large mixing angles, up to a « 1, are compatible with the data, but smaller 
mixing angles are slightly preferred. The strongest constraints on the mixing angle comes 
from measurements of rzz and ry/w since those rates decrease as cos 2 a with increasing a. 

So far we have kept cq = which is favored by the global fit. Obviously for very large 
mixing angles a ~ n/2 this leads to disagreement with the data, since r Z z,ww ~ > in this 
regime. To estimate the viability of the pure CP-odd scenario, we can instead fix a = n/2 
and redo the fit. Since T(0 — > bb) is zero here, the dominant decay of is into gluon and 
photon pairs and into Zj, while the WW and ZZ decays are suppressed by the three body 
phase space. A realistic value for both r 77 and rzz then requires an approximate cancellation 
of the effective photon coupling, (cqCb + s 2 d cw) 2 <C c 2 B , c^y. 

It is still impossible, however, to achieve rzz ~ 1 if the operator coefficients are restricted 
to |cj| < 47T, as suggested by perturbativity. While the measured r 77 is well reproduced, the 
model is excluded at the 99.7% level due to the absence of a signal in all other channels. 
Note that allowing for nonzero x u , Xd, i-e. allowing the CP-odd component to couple to 
SM fermions, does not lead to an improved fit. The reason is that while a nonzero x u can 
increase the production cross section, it also increases T{h — > gg) and therefore the total 
width of so that the effects drop out in the ratios rzz and ry/w- We therefore arrive at 
the very strong conclusion that a pure CP-odd resonance is excluded at the 3a level in any 
perturbative extension of the SM. 

Finally note that our setup also allows us to study the fermiophobic limit x u ^ = y u ,d = 0. 
In the absence of mixing, this parameter point disagrees with the data at the 3a level. Once 
mixing and nonzero coefficients for the dimension five operators are allowed, we instead find 
a good fit to the data for 

a = 0.84, c G = 0.94, c B = 0A7 , c w = 0.16 , (27) 

with a x 2 similar to the SM fit. The Tevatron evidence for Higgs to bb decays is in tension 
with this parameter region, however stronger evidence of Higgs decays to SM fermions is 
required to probe a mixed fermiophobic scalar. 

5 Future projections 

It is therefore worth estimating how much the constraints on CP violation in Higgs mixing 
can be improved with future data from the LHC, using only the rate information. We expect 
that a full analysis of the 2012 8 TeV data set will lead to a reduction of the uncertainties 
by roughly a factor two. 

In the long term, measurements in the 14 TeV run will not only improve the sensitivity 
in the current search channels, but will further add rate measurements in not yet observed 
channels. Specifically we have considered the following additional channels from Ref. [27] : 

• h-> WW* (VBF tag) 
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Figure 5: Projected sensitivity on the CP mixing angle a and CP-odd top Yukawa coupling 
x n from upcoming LHC data on Higgs rate measurements. The plot is similar to Fig. [IJ 
under the assumption that all rate measurements have a central value consistent with the 
SM. Shown are the 95% C.L. limits for current errors but SM-like central values (light solid), 
quadrupled statistics per experiment at the end of the 8 TeV run (light dashed), and expected 
errors for 300 fb~ x at 14 TeV, including only the channels discussed in L 271 and combining 
two experiments (green, thick). 

• h — > /jl/jl (inclusive) 

• Vh, h — > 77 

• tih, (h —> 77 and h —> yU/i) 

In addition the note specifies projected sensitivities for the ZZ, WW (inclusive), 77 (inclu- 
sive), 77 (VBF) and rr (VBF) channels, which are also included in our estimate. Projections 
are not given for the bb and Z7 channels. 

Of course it is impossible to predict if the central values will remain the same or shift 
when more data will be analyzed. For concreteness, it is interesting to estimate how well 
the mixing angle can be constrained under the assumption that the rate measurements will 
converge towards the SM predictions. To illustrate this, in Fig. [5] we show three curves in 
the a-x u plane. As in Fig. CD we set Xd = 0, y u = yd = 1, and cq = cw = cb = 0. However 
instead of using the currently measured signal strengths, we assume that rx,measured — 1 f° r 
all channels. The three curves then are the 95% C.L. limits given (i) the current error bars 
(light solid), (ii) assuming a factor two improvement on the error bars (light dashed), and 
(in) expected uncertainties based on the 300 fb _1 14 TeV ATLAS projections and including 
an additional factor of two in the statistics assuming that CMS reaches similar sensitivity 
(green, thick). Note that for the last curve, we only include the channels for which projections 
are given in [27J. 
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As evident from the plot, the mixing angle a in this scenario could be constrained to 
a < 1.1 at the end of the 8 TeV run and to a < 0.7 with data taken at full energy and 
luminosity. The bounds thus will become significantly stronger, although they are weakened 
by the dependence on other parameters, such as x u . 

With larger available data sets it becomes feasible to constrain CP properties from an- 
gular distributions, and the projected sensitivity at the end of the 8 TeV run is similar to 
Fig. [8]. There are currently no estimates for the angular analysis of decay products of a 
125-GeV boson at 14 TeV, but it is likely that it will lead to superior limits compared to 
Fig. [51 since this method is not affected by the a priori unknown couplings x u ^. 

6 Conclusions 

The CP properties of the newly discovered boson with mass in the range 125-126 GeV can 
already be constrained with existing data on production rates and branching ratios. The 
main reason is that a CP-odd pseudoscalar generally has suppressed couplings to W and 
Z bosons, since they are generated only by higher-dimensional operators. In this paper, 
the 125-GeV resonance has been assumed to be a general mixture between a CP-even 
and CP-odd scalar field, and bounds on the mixing angle have been derived in a variety of 
different scenarios: 

(i) fermion Yukawa couplings fixed to their Standard Model values; 

(ii) the overall scale of up- type and down- type Yukawa couplings to the CP-even and 
CP-odd components may float freely and independently; 

(Hi) the signal peak near 125 GeV is comprised of two particles, and 0', which are the 
two mass eigenstates of the mixed CP-even and CP-odd scalar fields; 

(iv) addition of higher-dimensional operators that couple the CP-odd component to gauge 
bosons; 

(v) a special case of (iv) with vanishing Yukawa couplings (fermiophobic scalar). 

Using the most recent Higgs search results released by ATLAS and CMS in July 2012 [H[2], 
it turns out that the possibility that is a pure pseudoscalar is already excluded at more 
than three standard deviations, assuming that the new-physics sector is weakly coupled. 
Nevertheless, large values of the CP-mixing angles, a ~ 1.0... 1.3, are still allowed that 
the 95% confidence level, although smaller values are in better agreement with the data. 
Interestingly, a non-zero mixing angle, a ~ 0.15... 0.4, together with modified Yukawa 
couplings, produces a slightly better fit than the Standard Model, although the difference is 
not significant. If one allows the possibility that two mixed scalars contribute to the observed 
resonance peak (scenario (Hi)) no conclusive constraint on the mixing angle can be derived 
from current data. 

Besides the mixing angle, meaningful limits on the Yukawa couplings and couplings 
strength of higher-dimensional operators of have been obtained in all scenarios. The bound 
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on these parameters and a are expected improve significantly with increased statistics, which 
will lead to more precise measurements of the production rates and branching ratios, as well 
as open up the possibility to constrain the CP properties by studying angular distributions 
of the decay products. 
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Appendix 

The Lagrangian (00) leads to the following Feynman rules for the coupling of the CP-odd 
scalar A to SM gauge bosons: 

^7 %t)t CW e ^ {P2 ' P ' )M - (28) 

AZZ S ^)i CW e " UPU{P2 - P3)M ° ' (29) 

AlZ SJ ^^^(P^PMPA, (30) 

AWW _^L_ e ^ (p2 _ , (31) 

AGG *L_ e ^ {p2 _ pj)^), , (32) 

where and Pi,2,3 ar e the momenta of the first, second and third particle flowing into the 
vertex. 
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